Varying the effective refractive index to measure optical transport in random media 
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We introduce a new approach for measuring both the effective medium and the transport prop- 
erties of light propagation in heterogeneous media. Our method utilizes the conceptual equivalence 
of frequency variation with a change in the effective index of refraction. Experimentally, we mea- 
sure intensity correlations via spectrally resolved refractive index tuning, controlling the latter via 
changes in the ambient pressure. Our experimental results perfectly match a generalized transport 
theory that incorporates the effective medium and predicts a precise value for the diffusion constant. 
Thus, we directly confirm the applicability of the effective medium concept in strongly scattering 
materials. 

PACS numbers: 42.25.Dd, 78.20.Ci, 78.67.-n, 42. 30. Ms. 
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Light propagation in multiple-scattering media is dom- 
inated by speckle, a highly irregular intensity pattern de- 
pendent upon spatial (or angular) or time (or frequency) 
coordinates and brought about by interference. Corre- 
lations, which are inherent properties of speckle despite 
the apparent irregularity, provide important information 
about transport parameters. In fact, the description of 
intensity correlations is at the heart of understanding 
transport theory 1]. In the past, correlations have been 
measured in time 0] and frequency 0] as a means of, for 
example, determining the diffusion constant of light. 

What is often not fully appreciated is the degree to 
which effective medium properties are essential for de- 
termining correlation functions, and thereby transport 
properties. When developing a theory of transport, the 
effective medium is an essential building block [J]. In 
practice, fundamental problems arise when one attempts 
to access the effective medium properties of strongly scat- 
tering materials. In a turbid medium, the coherent prop- 
agation, upon which the effective medium properties are 
defined, decays on the scale of a few mean free paths 
and becomes unmeasurable in the overwhelming bath of 
diffuse intensity. On the other hand, the assumptions 
needed for developing standard effective medium models 
do not hold for such materials, which makes the mod- 
els unapplicable. These difficulties raise the question of 
what exact role the effective medium plays in determin- 
ing transport properties in strongly scattering materials. 

In this letter we first highlight the central role of the 
effective refractive index in transport theory. We then 
propose and measure a new correlation function that di- 
rectly incorporates the effective medium approach. We 
validate this approach by showing the precise agreement 
between the theory and experimental results. This com- 
parison requires defining a new directly measurable effec- 
tive medium parameter, which we call the "tuning re- 
sponse", as 



The parameter 5 relates the change of the effective re- 
fractive index n e to the changes in the refractive index 
nh of one of the components, which can be easily exper- 
imentally controlled. We measure the precise value of 5 
with our method and use it to extract the diffusion con- 
stant. Finally, we use 5 to test several standard models 
of the effective medium. 

Our measurement method is based on controlled 
change of the optical path length distribution by refrac- 
tive index tuning (RIT). It is the optical analogue of a 
class of experiments in condensed matter physics that 
control the elongation of electron trajectories by using 
magnetic field. Those experiments resulted in the obser- 
vation of electronic weak localization and universal con- 
ductance fluctuations The concept of changing the 
path length distribution has also been exploited for dif- 
fusing wave spectroscopy [g, M| , and for other types of 
waves in evolving media [il li^liol]. Given the predom- 
inance of disordered photonic media in a wide variety 
of fields including biology 11| , advanced materials [l2| , 
solar cells [13], and in general modern photonics, we ex- 
pect our approach to have broad cross-disciplinary appli- 
cation. 

In standard transport theory Q , the effective refractive 
index shows up in the early stages where the the averaged 
amplitude Green's function is introduced: 
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where t s is the scattering mean free path. This form of 
the Green's function takes care of internal resonances of 
scatterers but does not hold in a regime of strong spatial 
dispersion, for example in photonic crystals, or when the 
material is anisotropic. One can introduce the effective 
(complex- valued) wavenumber by K = 2^ + i-^jr- 

We introduce a new correlation function, which is cen- 
tral to our measurements, defined as: 
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C u , u+ Au{n e ,n e + Art e ) = N [{I u (s;n e )Iu+Au(s;n e + An e )) - (I u (s; n e )) (I w+ Au>(s; n e + An e )>] , 



(3) 



where I u (I) is the far-field specific intensity at direc- 
tion s and (.) denotes averaging over a narrow spread of 
directions. The normalization constant N is fixed by re- 
quiring Cu,^(n e ,n e ) — 1. Intensity correlations are typ- 
ically classified into different contributions: C^ 1 - 1 , C^ 2 \ 
and Their relative magnitude depends on the 

scattering strength of the material. We will present the 
closed form for the first order contribution G^ 1 ' , which 
is sufficient to describe our measurements, but our argu- 
ment holds for the higher-order terms as well. 

The calculation of correlation function for the case 
of Arte = can be found in many papers and text- 
books P, [3] • The actual expressions depend on the ge- 
ometry. Closed forms have been presented for the case of 
an infinite medium [l5| , semi- infinite medium [l6| , and a 
slab [3| . Yet conventional definitions of this type of cor- 
relation function have not included variations of n e . The 
crucial point of the theoretical part of the present paper 
is that the formulae calculated for frequency correlations 
can easily be generalized to the case of An e ^ 0. 

In our regime of consideration where 
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we can neglect variations of ImK in the average Green's 
functions. As a result, the change in the average Green's 
function © is only dependent on the product n e uj. That 
is to say changing the frequency by Alu is equivalent to 
changing n e by An e if n e Au = LoAn e . 

We prove that the same symmetry between frequency 
and effective refractive index variation holds for the in- 
tensity correlation function ([3]). The basic ingredients for 
deriving this correlation function are the average Green's 
function |(3J) and an irreducible scattering vertex We 
symbolically denote the irreducible scattering vertex as: 
U(u),u> + Alu; n e ,n e + An e ), and assume that it is local, 
thus we have dropped the three momentum variables. 
Changes caused by U are compensated by those result- 
ing from AlmK due to energy conservation. This cance- 
lation holds even without considering condition (TJ|). As 
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FIG. 1: The experimental setup. The sample is placed in 
pressurized chamber and illuminated by a HeNe laser. A sin- 
gle polarization of the transmitted speckle is recorded by a 
CCD-camera. 
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FIG. 2: RIT autocorrelation coefficient as a function of Ann- 
Data points show measurement results for 3 slabs of polyethy- 
lene filters (L = 1, 1.5, and 2 mm for black squares, red circles, 
and blue triangles, consequently). The lines show the fit to 
theory (Eq. ©). 



a result, replacing n e (uj + Alu) by (n e + An e )(uj + Alu) 
does not affect any algebraic step in the derivation of 
Cuj.uj+Auj(n e , n e ) [141 and thus yields the generalized cor- 
relation function C W)W +A w (H e , n e + An e ). 

To take one example consider a slab of porous ma- 
terial, for which the host refractive index, rih, may be 
tuned. The correlation is then measured as a function 



of rih instead of n e 
given by: 



In this case, Cu,u(nhi n h + An/,.) is 



C£l(n h ,n h + An h ) = 



T S An h 



cosh \JrsArih — cos \Jts Ann 



(5) 



with RIT decay coefficient (generalized analogue of dif- 
fuse decay time) 



Tg 
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where we have used the tuning response S defined in 
Eq. ©. 

We apply our method to typical samples, which are 
representative for a a large class of multiple scattering 
materials. The sample dimensions are selected such that 
they can be characterized by both RIT and time-resolved 
measurements, allowing us to show the accuracy of our 
technique. We have used a composite material with open 
channels consisting of a solid backbone and a gaseous 
host. The refractive indices of the solid backbone and the 
gaseous host medium are n s and respectively. In our 
measurements the index of the host is tuned by increasing 
the pressure of the gas. 

As samples we used slabs of commercially available 
porous plastic air-filters (XS-7744, Porex Corp.) avail- 
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FIG. 3: Measured relative change of n e as a function of Ann- 
The value of An e /n e is extracted by measuring this spectral 
shift. The refractive index of air is calculated based on the 
Edlen formula [23|. Inset: The shifted cross-correlation co- 
efficient C versus the spectral shift Av for certain pressures 
(From right to left: AP = 0.0, 0.075, 0.18, 0.28, 0.37, and 
0.46 bar). The peak position of each curve denotes the spec- 
tral shift, which is proportional to the relative change in the 
effective refractive index. The spiky feature at zero position 
is an artifact of the slight inhomogeneity over the detection 
efficiency of CCD pixels. 



able in different thicknesses. They are composed of sin- 
tered polyethelene spheres with a broad size distribution 
of 7 to 12 /im. The refractive index of the polyethelene is 
n s — 1.49. This material has a porosity of <f> = 0.46±0.02 
and a mean free path of t a — 20.6 ± 0.2 /im in the fre- 
quency range of interest. The porosity was determined 
by weighing a larger sample of the same material. The 
mean free path was determined from total transmission 
measurements. These measurement also revealed a weak 
dependence of £ s on frequency, justifying the assump- 
tion (TJ). 

The sample was kept in a pressurized chamber and 
illuminated by a HeNe laser at 632.8 nm. Part of the 
transmitted speckle pattern was filtered by a linear po- 
larizer and recorded on a 16-bit CCD-camera with 10 6 
pixels. The recorded image consisted of ~ 10 4 indepen- 
dent coherence areas. Gradually tuning the air pressure 
in the chamber changes n^, causing the speckle pattern 
to evolve. This evolution allowed us to directly measure 
the autocorrelation coefficient defined in Eq. ^ . 

This experiment was performed on three samples. The 
extracted autocorrelation coefficients are plotted in Fig. [2] 
versus An/,. Each data series is fitted to the correla- 
tion function of Eq. ([5]) with a separate single fitting pa- 
rameter, tj. We see an excellent agreement between the 
theory and the experimental results. This agreement is 
an important confirmation of our theory which is built 
upon the concept of effective wavenumber in random me- 
dia. This concept has never been examined so directly in 
strongly scattering media, simply because it was a very 
difficult parameter to access experimentally. 

Our experimental results not only match the theo- 
retically predicted functionality, but are quantitatively 



precise. We show this precision by extracting the dif- 
fusion constant from the RIT decay coefficient t$ and 
comparing it with the diffusion constant measured by 
an standard method. This requires the value of 6 
which we can measure using the symmetry relation dis- 
cussed earlier in the paper. This symmetry implies that 
Cuj,uj+Aui(n e ,n e + Art e ) has a peak equal to one for a 
nonzero shift Alo given by 

[n e + An e ) ( + Au> ) = n e u> or = . (7) 

to n e 

In other words, S can be extracted from the spectral shift 
of every individual speckle spot while the host index of 
refraction rih is varied. 

To monitor this spectral shift the experimental setup 
was changed. The light source was replaced by a 
white-light super-continuum laser (Fianium). The CCD- 
camera was replaced by a spectrometer, which was run in 
the imaging mode. The entrance slit of the spectrometer 
selects a transmission direction in form of a narrow rect- 
angle, which contains roughly 10 independent coherence 
areas along the slit. The beam is spectrally resolved per- 
pendicular to the slit direction by a grating. This config- 
uration allows us to simultaneously monitor the spectral 
evolution of several speckle spots while rih is changed. 

The measurement was performed on a 100-/im-thick 
slice of the polyethelene filter. A thin sample was nec- 
essary due to the limited resolution of our spectrometer. 
The intensity spectrum between 649 and 651 nm was 
measured while changing the pressure in the chamber 
from 1 to 1.5 bar in steps of 10 mbar. 

RIT correlation functions at different frequencies are 
calculated from the measurements and plotted as a func- 
tion of Auo for each pressure. A collection of these 
plots is shown for six different pressures in the inset of 
Fig. [3] as typical representatives. As predicted, the peak 
value is shifted and equals unity within the experimen- 
tal error. Using our symmetry relation ([7]), the relative 
change of effective refractive index is calculated and plot- 
ted as a function of An„ in Fig. We get a value of 
5 = 0.212 ± 0.003 from the slope of a linear fit to the 
data, for these set of samples, with a remarkably high 
precision. 

In principle, it is possible to extract the diffusion con- 
stant from each separate RIT measurement. To im- 
prove the accuracy, we plot the auxiliary parameter 
a = AitvS/ts versus L~ 2 from the three measurements 
presented in Fig. [5J This plot is shown in Fig. 2J The 
slope gives the diffusion constant D = 2.2 ± 0.1 m 2 /ms 
using Eq. ([5]). 

Next we compare our result with the diffusion con- 
stant measured by conventional time-resolved method. 
The diffuse transmission through the same samples was 
recorded using a time-correlated photon counting setup 
and a subpicosecond pulsed laser at 600 nm wavelength. 
The resulting temporal decay curves are presented in 
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FIG. 4: The measurement parameters a 



and b : 



are plotted versus the inverse of the thickness squared. Black 
squares and red dots correspond to the index tuning and time- 
resolved measurements respectively. The slope of each data 
set is equal to the diffusion constant D measured by each 
specific method. The fact that the two slopes are equal proves 
the consistency of our index tuning method with the time 
resolved measurements. 



the supplementary information. Following the guidelines 
of [l^| and considering the response function of the detec- 
tor, the diffuse decay time Td = L 2 /ir 2 D is extracted from 
a fit to these measurements. These results confirm that 
the absorption has a negligible effect on our experiments. 
The diffusion constant is calculated from the slope of a 
linear fit of the parameter b = -^r- vs. L~ 2 (Fig. Q|, 
D = 2.1 ± 0.1 m 2 /ms, which is in excellent agreement 
with the result from the index tuning experiment. 

Our precise measurement method for measuring 5 is 
applicable to a variety of samples, and assumes no specific 
effective medium model. The parameter 8 can also be 
predicted using models of the effective medium 
22] which are all based on long- wavelength limits. 

For illustrative purposes, we present the predicted 
value for n e and S by three of the more popular effec- 
tive medium theories: (i) the average permittivity 24J, 
(ii) the Maxwell-Garnett and (iii) the Bruggeman mod- 
els. The predictions are n e = 1.29 ± 0.02, 1.25 ± 0.02, 
and 1.22 ± 0.02 and 6 = 0.28 ± 0.02, 0.45 ± 0.03, and 
0.52 ± 0.03 respectively. These predictions for S differ 
from each other and from our experimental result by as 
much as a factor of two. It is also worth noting that the 
predicted value of S is more model-dependent than n e 
itself. The main reason for this discrepancy is perhaps 
the assumption of the long wavelength limit, which is not 
reached in our experiments. Our unambiguous measure- 
ment of S highlights the necessity of more sophisticated 
models of the effective medium [2l| . 

In conclusion, we have presented a new effective 
medium quantity and a new transport correlation func- 
tion and shown how to precisely measure both of them 
by using refractive index tuning. Our measurements di- 
rectly test and approve the validity of assuming an effec- 
tive wavenumber for an inhomogeneous medium, which is 



very important for describing all sorts of photonic meta- 
materials. Using these two quantities one can measure 
several dynamic transport properties with high precision, 
as we have demonstrated for the lowest order cor- 
relations and the diffusion constant. Measuring higher 
order correlations and studying Anderson localization is 
a natural follow-up to our research. 
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n e = (f>n h + (1 - <f>)n s 
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FIG. 5: Time resolved intensity decay curve for slabs of polyethylene niters. Different datasets correspond to different sample 
thicknesses (L = 1, 1.5, and 2 mm for black squares, red circles, and blue triangles, respectively). The lines show the fit to the 
diffusion model after considering the absorption and the response of the detector. 
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